Report for Model 6

— for both jointly and separately

2024-08-19

1. Joint Modeling

1.1 Trace plots for convergence check

The current MCMC setting is:

  • 300,000 iteration;
  • 200,000 burn-in;
  • 100 thinning.

1.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum      1.002       1.01
## LDevsum      0.999       1.00
## dh0          1.918       3.57
## dh1          1.668       3.63
## dh2          1.135       1.45
## dl0          1.658       2.87
## dl1          1.334       2.03
## dl2          1.074       1.29
## dl3          1.258       1.85
## 
## Multivariate psrf
## 
## 1.49

1.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

1.4 WAIC results

LevelH LevelL
DIC 1231.31073 22756.8074
DIC3 1164.66422 20662.5287
PWAIC 46.58249 994.8117
WAIC 1194.22015 21039.8413

2. Separate Modeling

2.1 Trace plots for convergence check

The current MCMC setting is:

  • 300,000 iteration;
  • 200,000 burn-in;
  • 100 thinning.

2.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.00       1.01
## LDevsum       1.00       1.02
## dh0           1.23       1.79
## dh1           1.12       1.40
## dh2           1.02       1.09
## dl0           1.11       1.12
## dl1           1.15       1.51
## dl2           1.03       1.11
## dl3           1.19       1.37
## 
## Multivariate psrf
## 
## 1.21

2.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

2.4 WAIC results

LevelH LevelL
DIC 1377.52079 23110.679
DIC3 1247.01158 20729.655
PWAIC 90.84391 1040.447
WAIC 1316.35844 21133.344